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Based on the Keldysh formalism, the Boltzmann kinetic equation and the drift diffusion equation 
have been derived for studying spin polarization flow and spin accumulation under effect of the time 
dependent Rashba spin-orbit interaction in a semiconductor quantum well. The time dependent 
Rashba interaction is provided by time dependent electric gates of appropriate shapes. Several 
examples of spin manipulation by gates have been considered. Mechanisms and conditions for 
obtaining the stationary spin density and the induced rectified DC spin current are studied. 
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I. INTRODUCTION 

Spin transport in semiconductor heterostructures has 
recently attracted much attention owing to a perspective 
of its practical application for quantum computing and 
communicationsiSi^ On the way of fundamental studies 
and applications of the spin transport the key problem 
is how to generate, detect, and manipulate the electron 
spin polarization. Many ideas have been proposed to 
achieve control of the spin using magnetic materials, ex- 
ternal magnetic fields and optical excitation (for a review 
see Ref. |3) ■ At the same time, a challenging goal retains 
to employ only the electric control integrated by means 
of gates into a high mobility semiconductor heterostruc- 
ture. 

A promising opportunity on this way opens due to 
the strong Rashba spin-orbit interaction (SOI) in nar- 
row gap two-dimensional (2D) electron systems.^ An im- 
portant feature of this interaction is its tunability which 
can be achieved by varying the gate voltagei^ The field 
effect transistor was the first proposal utilizing this phe- 
nomenon J. The gate control of the spin current employ- 
ing the Aronov-Casher effect^ was considered in Ref. 0. 
The electric dipole spin resonance controlled by the time 
dependent gate was studied in Ref. ^3 

The SO interaction effect alone, or in combination with 
the external electric field, allows to approach an impor- 
tant goal to create spin currents and spin polarization by 
entirely electrical means, not involving any of the mag- 
netic materials or optical excitation. One of the exam- 
ples is the recently predictedii (see also Ref. IT^ and 
observediSii^ spin-Hall effect in 2D and 3D electron and 
hole gases, where the spin current is driven by the elec- 
tric field. The closely related phenomenon is the spin 
polarization of 2DEG in response to the parallel electric 
field. "'^^ Another method utilizes the time dependent gate 
to modulate the shape of a quantum doti& or the strength 
of the SO coupling constant in IQiLiS and 2D systemsii^ 
In the latter case an efficient spin current generation can 
be attained in the presence of high frequency (hundreds 
MHz or higher) gate-voltage variations. The physics of 



this phenomenon is simple. The Rashba spin-orbit in- 
teraction resembles an interaction with a spin dependent 
gauge field. Therefore, its time derivative produces a mo- 
tive force on electrons, similar to the electromotive force 
from the time dependent electric vector potential. The 
important difference is, however, that the force created 
by SOI acts in opposite directions for oppositely polar- 
ized spins. Such a method of spin current generation is 
convenient for implementation in conventional semicon- 
ductor heterostructures. By a proper choice of the gate 
shape it allows to create and rectify the AC spin current 
or accumulate the spin polarization at a given location. 
However, the theory in Ref. is restricted to the spa- 
tially homogeneous case. Therefore, within this theory 
one can not properly consider boundary effects, as well 
as the spin current generation due to a time dependent 
gate of the small area. Also, it is impossible to study 
any effects of the spin current generation accompanied 
by spin accumulation. 

In the present work we develop a theory which is based 
on the Boltzmann transport equations. This approach is 
quite universal and it allows us to study the spin trans- 
port under an arbitrary space-time dependent SO inter- 
action, providing that a characteristic scale of this depen- 
dence is within the range of applicability of the semiclas- 
sical approximation. This means that the spatial varia- 
tions of the SO interaction are assumed to have the scale 
Ar larger than the electron Fermi wavelength and the 
scale of its time variations At ^ h/Ep, where Ep is the 
Fermi energy. For even more larger time-space scales, 
such as Ar I and At » r, where I and r are, respec- 
tively, the mean free path and the mean free scattering 
time, we use the Boltzmann equation to derive the drift- 
diffusion equation. Within this theory we consider the 
following problems: 

1. Spin current generation by a finite size time de- 
pendent gate. In this case the spin polarization is 
pumped by the gate into 2D regions adjacent to 
it. This polarization further diffuses away from the 
gate, as well as in the backward direction. 
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2. Spin current generation in a gas confined within an 
infinite 2D strip. In such a geometry the generated 
spin current with the polarization perpendicular to 
boundaries flows freely along the strip. At the same 
time, spins polarized parallel to the boundaries are 
accumulated near them, opposite spins near oppo- 
site banks. 

3. Rectification of the AC spin current. It will be 
shown that the spin polarization generated by a fi- 
nite size AC gate in some cases contains the DC 
component due to periodic variations of the elec- 
tron density under the gate. The DC component 
appears as a result of interplay of two periodic pro- 
cesses: oscillations of the SO coupling constant and 
oscillations of the 2D electron density. 

This paper is organized as follows. In Sec. II, the 
Boltzmann equation is derived for calculation of the spin 
transport in the presence of the time-dependent SO in- 
teraction. In Sec. Ill, we apply the Boltzmann equation 
to some of the spin transport problems. The derivation 
of the drift-diffusion equation and its applications are 
presented in Sec. IV. The results of this work are sum- 
marized in Sec. V. 



II. BOLTZMANN TRANSPORT THEORY 



e±(k) 



2m* 



± aok. 
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For 2DEG density of interest the Rashba spin-splitting 
energy Aq — 2aok will be assumed to be small compared 
to the Fermi energy, Aq <C Ep. In the following, we set 
h = I for convenience. 

Let us consider the time-modulated nonequilibrium 
electron system at temperature regime T <^ Ep, assum- 
ing the semiclassical conditions to be satisfied. Hence, 
the time intervals At and distances Ar over which 
all quantities vary significantly satisfy the inequalities 
EpAt ^ 1 and ppAr ^ 1. These conditions are nec- 
essary for derivation of the kinetic equation. We start 
this derivation from introducing the matrix of nonequi- 
librium Green's function in the Keldysh spacei^S 



G-sil,2) G-+(l,2) 



G:^"(1-2) G:+(1,2) 



(3) 



where 1 and 2 denote two points ri, ti and r2, i2, in the 
time-space and a, /? are spinor indices. We will assume 
that besides the time dependent forces in the Hamilto- 
nian (Q, the electrons are also subject to scattering from 
randomly distributed static impurities. This scattering 
can be taken into account by introducing the self energy 



We consider spin transport in a noninteracting 2DEG 
confined within a narrow-gap semiconductor quantum 
well (QW). The 2DEG is applied atop by a gate under 
time dependent bias. The inversion asymmetry of the 
quantum well confining potential is controllable by an 
external gate, and hence electrons experience a tunable 
Rashba-type SO interaction.-^ The system is therefore 
described by the effective mass Hamiltonian 



H{t) = 



2 m 
+U{v,t), 



M-i[a(r,t)p-f pa(r,t)] 



(1) 



where m* represents the effective mass, p — {px,Py) is 
the electron momentum in the 2DEG plane, and a(r, t) 
denotes the time-dependent Rashba coupling constant. 
a(r, t) contains both a static part ao and a dynamic part 
ai(r,t) due to the time dependent gate. In addition to 
time dependence, ai can also vary in space depending 
on the gate configuration and applied bias. We denote 
M — z X (T where z is the unit vector along the growth 
direction and cr = {cFx, cry, a^} is the vector of the Pauli 
matrices. Further, U (r, t) indicates the external poten- 
tial energy of an electron in the electric field produced 
by the gate. The static part of the Rashba coupling con- 
stant lifts the spin degeneracy resulting in a momentum- 
dependent spin-splitting of the conduction band, namely 
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y+- y++ 
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The corresponding matrix Dyson equation is of the form 

dM3 [^°^(l,4)i;^5(4,3)^i;3(3,2)] (5) 

or the conjugate form 

^0/3(1,2) = ^°M1'2) 

d3dA [^„^(l,3)i;^^(3, 4)^0^(4, 2)] ,(6) 

where the functions in the integrand are matrices both in 
real space and in spin space being combined by the rule 
of matrix multiplication. Acting on the Dyson equation 
from the left (or from the right on its conjugate) by the 
operator 



d 



(7) 



where the suffix j = 1 or 2 indicates the differentiation 
is with respect to the variables tj and Vj. After some 
algebra, we obtain the equation^S 
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a(3 



a(2)V2 • (5^7(1, 2)M^0) + Vi • (m„^^^/3(1, 2) 
V2a(2) • (^a^M^^) + Via(l) • (m„^^^/3(1, 2)) 



(8) 
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which is equivalent to the set of integro-difFerential equa- 
tions for the component Green functions. Here the scat- 
tering integral is defined bjiSS 



(isc) 



Q/3 



= / d3 



fi;„^(i, 3)^^^(3,2) 



-^„^(i,3)i;^/3(3,2)7 



(9) 



where f is the matrix, such that r++ = — t = 1 and 
T+- =T-+ = 0. 

In the semiclassical regime, it is convenient to trans- 
form the variables to the Wigner coordinates: r = ri — r2, 
R = i(ri -t- rs); t = h - t2, T ^ ^{h + ^2)- The differ- 
ence variables r, t vary on a microscopic scale, while the 
center variables R and T are macroscopic variables. In 



the quasiclassic approximation the Green functions and 
self energies vary slowly with respect to the center vari- 
ables. Accordingly, only linear gradient expansions will 
be taken into account in JSJii^It is appropriate to define 
the space-time Fourier transform to the fast variables 



X exp [i (k • r — ojt)] 



(10) 



where the uppercase indices i,j take the values "4-" or 
"— " . Denoting as G*-' (R, T; k, oj) a matrix with elements 
G'^^(R, T; k, w) and making the gradient expansion of 
© one can write the quantum kinetic equation in the 
form 



dT 



V • VrG^^ 



[k • M, G'- 



-a 



{M,VrG^^} 



-ivKa.{k.M,VkG'^} + i|^|k.M, 



dG'^ 



- VU ■ VkG'J 



dU dG'i 



(11) 



where v = k/m* and {A, B} denotes the anticommuta- 
tor. The scattering part of this equation is obtained from 
(|^ in the leading order of the gradient expansion with 
respect to the "slow" R, T variableSfSP In calculations be- 
low we will employ the three Green functions. There are 
G ^, retarded G^ and advanced G° functions. G*" and 
G'^ are defined by the equations 

r^r — h r^+— /^++ 

G" = G-+ ~G++ ^G-- ~G+-. (12) 

The same equations are valid also for the retarded and 
advanced self energies with the only difference that 
and S ^ enter with the signs opposite to signs of G"' 
and G ^ in H12I) . In terms of these functions the corre- 
sponding scattering parts of (|11|) are written as 

I-+ = -E''G-+ + G-+S° - G'^E-+ + S-+G^ (13) 



/,'; = -[S^G'■];/^, = -[E^G^ . (W) 



In order to determine the self energy we assume, 
for simplicity, that the impurity scattering potential is 
isotropic, spin independent and short range. Hence, it 
can be simply written as V{r) = VS{r) so that the corre- 
sponding Born scattering amplitude does not depend on 
the electron wavevector. In the quasiclassic approxima- 
tion ignoring weak localization effects the averaged over 
random impurity positions self energy is thus given bySS 

^ (_i)i+J|y|2^G'^J(R,T;k,cj). (15) 

k 

One can see that the so defined self energy does not de- 
pend on the wavevector. In the thermodynamic equi- 
librium the unperturbed retarded and advanced Green 
functions are easily found from (|15|) and 

GS = GS^ = (co-i/o + *r)-\ (16) 

where Hq = {k'^ /2m*) + apM • k. The elastic scattering 
rate F is obtained from H15|l . Substituting (|15|l into H15|) 
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one easily finds F = ■n\V\'^NF^ where Np = m*/(27r) is 
the one particle 2D density of statesi^ The +" unper- 
turbed Green function can be found from the equationSfi 

G-+ = -np{u){G\, - Gil (17) 

where nF('^) is the equilibrium Fermi distribution. 

The 2x2 matrix Wigner distribution function is de- 
fined as —iG '"(R, T; k, cj). At the same time, the den- 
sity matrix which is used for calculations of observable 
physical quantities, for example, spin density and spin 
current, can be obtained from the Wigner function by 
integration over w. Usually, the Boltzmann equation for 
this function is obtained by integration over uj of the 
quantum kinetic equation, such as the " — h" component 
of Eq. idJ. In our case, however, this method can not be 
used because of the term proportional to da/dT, which 
vanishes after integration. Moreover, it becomes a not 
simple task to write the scattering part in terms of the 
density function. Similar problem arises when one de- 
rives the Boltzmann equation for the system driven out 
of equilibrium by the electric field represented by a time 



dependent vector potential. This difficulty is resolved 
by shifting variables from wavectors to kinematic mo- 
menta or to velocitiesi^ In the problem considered here 
such a shift is not much helpful, because the velocity 
k/rn,* -|- Q!(R, T)M is a matrix in the spin space, not a 
number. We will employ a different method. Let us rep- 
resent the Wigner function in the form 

-iG-+ ^ inF{uj){G'' - G'') + F. (18) 

At C/(R,T)=0 and a(R,r) = ao the first term in the 
right-hand side turns into the unperturbed Wigner func- 
tion. Hence, the _F-function is not zero only due to devi- 
ation of the system from the original homogeneous ther- 
modynamically equilibrium state. Assuming this devia- 
tion to be small, we will linearize Eq. pi|l . omitting all 
products of F with terms containing the time and space 
derivatives of a and U . After substitution of Eq. H18() 
to the " — h" component of Eq. Ullfl , taking into account 
Eq. H13(l and the corresponding equations for G"" and G° 
with the scattering terms given by Eq. H14|l , we arrive to 
the following equation for F: 



OF 
df 

i da 



ia[k- M, F] + -a {M, VrF} 



G"} 



dnioj) 
du! 



^^(G^ 



G" 



dn{Lu) 



du! 



(19) 



where St[F] is given by 

St[F] = - FS" - G'^{F) + E(F)G^ (20) 

with E(F) defined as S(F) = \V\^ 'Zk F. 

A good approximation to the functions G^ and G" are 
the local equilibrium functions G[ and Gf defined by 

G^; = Gf ^ (u; - Hi + iT)-\ (21) 

where Hi = {k^/2m*) + a(R, r)M • k + U(R, T) - Ep- 
The local functions has the same form as the equilib- 
rium functions H16|l with the electron energy and spin 
splitting parametrically dependent on time an space. It 
can be seen from (^TJ, using and that 

the corrections to the local functions are proportional to 
the Green function gradients times the small parameter 
a/vp, where vp is the Fermi velocity. Such a small pa- 
rameter can be important for the spin-Hall effecliii, but 
not for the considered here case of the spin current driven 
by the time dependent SO coupling constant. One can 
further simplify Eq. Ijiyi) omitting proportional to this 
parameter terms in (|19|l . such as a{M, VrF}, which 
enters together with the much bigger v • VrF. After 
this simplification the linear operator acting upon F in 
the right-hand side of Eq. H19|) decouples into scalar and 



spin dependent parts. Further, since we are interested 
in the spin transport driven by the term in (|19|l pro- 
portional to the time derivative of a, one can ignore the 
term containing dU /dT, because it contributes additively 
into the linear response and gives rise to the spin-Hall ef- 
fect which will not be considered in the present work. It 
should be noted that the linearization is undertaken only 
with respect to the time derivatives, while local values of 
a and U entering into G\^°'^ can vary noticeably within 
macroscopic time-space scales. 

Now let us take a look at the scattering part H2U|) . One 
can easily see that since SOI is an odd function with re- 
spect to k, the self-energies S''^" calculated from (|15|) with 
the local equilibrium functions (|21|) are contributed only 
by the scalar parts of these functions. Therefore, omit- 
ting the small corrections of the order of A/Ep, where 
Ep is the Fermi energy, from (|21|l . H15|l and (|12|l we get 

^rp _ p^a ^ _2irF. (22) 

The other terms of (|20f) . the ones which contain S(F), 
can be also simplified. We notice that the spin depen- 
dent parts of G^^"'' contribute effectively with the small 
parameter A/Ep. Although small, these terms pro- 
vide coupling of the particle density represented by the 
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scalar Tr[S(F)] to the spin distribution function associ- 
ated with Tr[(TF] and vice versa. They are important in 
the spin-Hall efFec*SiS4. In our case we can ignore them. 
Taking into account all these simplifications we are ready 
to derive the Boltzmann equation for the particle distri- 
bution function in the space of particle coordinates and 
momenta. This function is obtained by integration over 
Lj of Eq. (O^" with C'" given by In the leading 

approximation we thus get 



I I — = riFZ + aM-k- 



27r 



dE 



+/k(R,T), 



(23) 



where npi = np {E + U (R.,T)) is the local Fermi dis- 
tribution function, where the coordinate and time de- 
pendent potential energy is added to the electron kinetic 
energy E. The function /k(R, T) is defined by 



/k(R,T)= j gi^(R,T;k,c.). 



(24) 



After integration Eq. H19|l over frequency we get the 
Boltzmann equation in the form 



St[/] - l^ + v VR/k + za(R,T)[k.M,/k] 



k M 



aa(R,T) dnpi 

df a^' 

The scattering term is obtained as 

k' 

where r = l/2r is the elastic scattering time. 



(25) 



(26) 



III. SIMPLE EXAMPLES OF SPIN 
TRANSPORT UNDER THE TIME DEPENDENT 
GATE 

In this section, we will employ Eq. H25|) to investigate 
electron transport properties of two-dimensional electron 
systems in the presence of time-dependent SOI. Two sim- 
ple examples of application of the Boltzmann equation 
will be considered: spin transport driven by a homoge- 
neous infinite gate and the ballistic transport due to a 
narrow time dependent gate. 

Before to proceed with these examples, let us define the 
spin current and spin density in terms of /k. According 
to definition of G '"^ the spin distribution in the space 
of particle coordinates and momenta is given by the spin 
dependent part of H23|l . Therefore, the spin density is 
obtained by integration of nonscalar terms in Eq. H23|) 
over k. Taking into account that the second term in 



Eq. H23|l turns to zero after averaging over k directions, 
we get the spin density 



PXR,r)-i^rr[fT7k(R,T)]. 



(27) 



The spin current definition is based on the one-particle 
spin-flux operator ■|{tT',v}, where the velocity operator 
V = V + (iM. Hence, using Eq. (|23|l . the spin current 
density can be written in the form 



dn 



Fl - 



dE 



4$]Tr [{v,a'}/k(R,r) 



(28) 



It is easy to see by a direct calculation that the first 
sum is zero. Therefore, only the nonequilibrium part of 
the spin current associated with /k(R, T) contributes to 
(|28|l . In this connection, it should be noted that the 
equilibrium spin current was found to be nonzero in a 
homogeneous 2D gasiSk This is due to quantum effects 
which are beyond the semiclassical approach used in the 
present work. 



A. Homogeneous case 

As a simple, yet nontrivial application, we consider 
the case of a large AC-biased gate such that the time- 
dependent region can be treated homogeneously. In this 
case one can omit Vn/k in Eq. H25|) and look for a so- 
lution of Eq. 1231 in the form /k = A(k)M. ■ k. Since 
/k = — /-k, the first term in Eq. H26|l turns to zero and 
the solution of l|25|l is easily obtained as 



/k = M • k 



a{n) 



dE '■ 



(29) 



where a(fl) is the Fourier transform of a(T) at the fre- 
quency O. For simplicity we have assumed that the elec- 
tron density, as well as dnpi/dE, do not change in time. 
On the other hand this dependence can be important for 
obtaining the rectified DC spin current. This effect will 
be discussed in the next section. Substituting (|29(l to 
(|28|l we find the spin current expression in accordance 
with Ref.lli 



1 pQ 



2n + 2iT 



(30) 



where p ~ 2^j.nF(A:) is the electron density, and the 
indices i and j denote the spin polarization and direction 
of the current flow, respectively. With the spin distri- 
bution function H29|l the spin polarization obtained from 
(|27|l is P = 0. Hence, in the homogeneous case no spin 
polarization is induced. 
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B. Current generation by a narrow gate 

We now consider the case of a narrow AC-biased gate, 
which is infinitely long in y direction while the width w 
in X direction is much smaller than the electron mean 
free path so that, except a small number of particles 
moving in the y direction the motion of electrons under 
the gate is ballistic. We are interested in the spin current 
flowing in the x direction. Since — —c^y, it follows 
from l|25|l that this current is polarized along the y-axis. 
Neglecting the scattering term we obtain from Eqs. H25|) 
and 



*2 



m 
47r 



dE 



da{x, T) dnpi 



dT dE 



2-K —iQ. + iqxV cos + 0+ ' 



(31) 



where qx is the wavenumber in the Fourier transform of 
the current with respect to the x coordinate and [..]n,ga; 
denotes the fi, qx Fourier component of the product in 
the square brackets. We assume that npi is coordinate 
an time dependent due to the gate electric potential. 
At Vlw <C w_F the frequency in the denominator of l|3HI 
can be neglected. Representing [iqxV cos + 0+)^^ as 
■n5{qx) / (y\cos(l)\) and returning to the coordinate and 
time representation we obtain 



Jn^,T) = -'^ J dx 



,da{x',T) 
'df 



kF{x,T), (32) 



where kpix' ,T) ^2fn*{EF - U{x',T)). The above 
expression is valid in the near vicinity of the gate, within 
the length of the electron mean free path. This ballis- 
tic result gives only a part of the spin current, the one 
associated with a direct generation action of the time 
dependent gate. It does not take into consideration the 
backflow of diffusion current due to the spin polarization 
accumulated on both sides of the gate. Such a diffusion 
current and the accumulated polarization will be calcu- 
lated within the drift diffusion theory in the next section. 
It is interesting to note that besides the AC component, 
the spin current also contains the DC component due to 
time dependence of fc^? in integrand of (|32|l . It can be 
easily seen that the DC current is obtained from (|32|l if 
harmonic oscillations of a(T) and fc_F(T) are phase shifted 
with respect to each other. Such a rectification effect will 
be studied in more detail in the next section. 



IV. THE DRIFT-DIFFUSION EQUATION 

In this section, we are interested in the time dependent 
spin dynamics in a disordered system, such that the char- 
acteristic frequency of the gate time dependence is much 
smaller then the elastic scattering rate 1/r, while the 
characteristic length of the spatial variation of /k(R, T) 



is larger than the mean- free path I. To this end, we start 
from Eq. I|25|) and represent /k in the form 



/k(R,T) = ^.gk(R,T)- 
Substituting ||221l into ^ we obtain 



^ + (V • VR)Sk 



2Q:ofc(Sk h) + hk 



-{Pe 

T 



9k) 



(33) 

da dnpi 
'df dE 

(34) 



where the unit vector h = (k x z)/fc and Pe — 
(27r/rn*) J2k, S{E' - E)g^,. With this definition of Pe, 
from H33|l and H27|) , the spin polarization can be expressed 
as 



m 
2^ 



dEPi 



(35) 



For simplicity we have assumed small variations of a and 
put a = ao in the third term of the left hand side of H34|l . 

In order to derive the drift diffusion equation the func- 
tion Qk is expressed from Eq. H34|l in the form 



9k = A 



- Ao(gk xh)-hk 



da dnpi 



dT dE 



(36) 



where A is the operator inverse to —iO. + v • Vr + 1/t. 
As we mentioned in the beginning of this section, in the 
diffusion approximation one can expand A^^ with respect 
to the small Vt and v • Vr compare to 1/r, so that 



A- 



T 1 + iVLt - t (v • Vr) + (v • Vr)' 



(37) 



The next step is to express g x h m the right hand side 
of (|36|l via Pe- It can be done by decomposing g into 
parallel and perpendicular to h parts, according to 

g = {gh)h + g^ (38) 

Taking the perpendicular projection of 1341) we find 



A{Pe xh)). (39) 



g^xh = -{A' + Ai)- 

T 

After inserting this expression into H36|l we use the ex- 
pansion H37|) and a similar expansion for (A^ -|- Aq)~^. 
It will be further assumed that Aq <C 1/t and only the 
terms not smaller than (Aqt)^ will be retained in the 
right hand side of Eq. (jSEJ. Since the current source 
in is proportional to dnpi/dE, for a degenerate 

electron gas the function Pe has a peak at the local 
Fermi energy Ep — C/(R, T). Hence, taking into ac- 
count the definition (|^ . this function can be represented 
as Pe = {2'K/m*)5{E -f t7(R,T) - Ef)P. Integrating 
(Eq. (|36|l 'l over energies and averaging over the angles of 
k we arrive to the diffusion equation for the spin trans- 
port 



dP 
'df 



(V • DV)P - 4Dm*ao [(z x V) x P] 
-r,P-r,P^z-i(Vxz)pr|^, (40) 
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where D = w|,;t/2 is the diffusion constant and Ts = 
D{rn*aQY is the spin relaxation rate. It should be noted 
that the diffusion constant can be coordinate and time 
dependent via the local Fermi velocity vpi. Except the 
last term on the right-hand side of Eq. (|40|l the above 
spin diffusion equation coincides with that derived earlier 
from the Green function formalism!^ In this equation the 
first three terms represent, respectively, spin diffusion, 
spin precession due to SOI and the D'yakonov-Perel— 
spin relaxation. The last term is the spin current source 
provided by the time dependent SOI. The spin relaxation 
gives the natural time scale Tjj = for the most of spin 
diffusion processes. One can also define the characteristic 
length of the spin density spatial variations as the spin 
diffusion length Ld — ^/DTd- The drift diffusion ap- 
proach is valid when Td ^ r. This condition is provided 
by the small Aq in comparison with l/r. 

For the following calculations we need an expression 
for the spin current. It can be found from its initial 
definition (|28|) . Inserting there /k written in the form 
^6'6\. the spin current is found as 



(41) 



Expressing according H36|) . from H37ll - 139|) we obtain 



1 ,,,-3 da 

2 dT 



(42) 



This equation shows that the spin current contains three 
components. The first and the second terms represent the 
usual diffusion current and the current associated with 
spin precession. These two contributions have been found 
in earlier worksi^ The third term is the new one. It 
represents the spin current generation due to the time 
dependent SOI. 

It is interesting to note that in the homogeneous case 
P = , and the spin current is simply of the form J'j = 
^e'^^pT{da/dT). This result coincides with Eq. (jSUJ) by 
taking the limit of il ^ in the denominator of Eq. (|30|l . 



A. The finite-size AC-biased gate 

We consider an AC-biased gate which is supposed to be 
infinite along y and with a width w wider then the mean 
free path I along a;, with a Rashba coupling constant 
a{x). The finiteness of the gate in the a;-direction results 
in X dependence of P. At the same time P does not 
depend on y. Retaining in Eq. (|40|l only derivatives with 
respect to a;, one can easily see that a pair of equations 
for P^ and P^ polarization components is decoupled from 
the equation for P^, and that the current source term 
enters only to the latter equation. Therefore, the solution 
of the diffusion equation is such that P^ = P^ = and 
the equation for the retained component is 

T.py + lsi.^T) (43) 



dT 



±D^py 

dx dx 



(a) 



(b) 



FIG. 1: Illustration of the spin current flow induced by a 
wide AC gate: (a) spin current as a function of x, which 
is uniform under the gate and exponentially decays outside 
the gate; (b) spin polarization as a function of x, which is 
accumulated at the gate edges. 



where the time-dependent spin-current source is given by 
S{x,T)^ prda/dT. 

We will apply Eq. H43I) to two problems. In the first 
problem the width of the gate will be assumed large, so 
that w ^ Ld- The opposite case oi w <C Ld will be 
considered in the second example. 



Let us assume that in the interval xi < x < X2 the 
diffusion constant is D(T) and a — aQ + a{T), while 
outside this interval D — Dq and a — ao. We will 
consider relatively slow variations in time of a{T), such 
that their frequency is much less than Fg. Consequently, 
the time derivative on the left-hand side of Eq. (|43|l will 
be neglected. The spin current induced in the range 
xi < X < X2 will be injected near edges at the points 
xi and X2, and the injected polarization will diffuse out- 
ward, as well as backward to the modulation gate region, 
as shown in Fig. 1. For analysis of this injection process 
it is enough to consider a vicinity of either edges. Let it 
be the right edge. Then, the solution of Eq. H43|l has the 
form 



py = Acxp 



\x - X2\ 

Ld 



(44) 



dx 



where Ld — (aom*)^^ . The coefficient A, in its turn, 
can be found from the spin current conservation at the 
boundary X2- On the right of the boundary the current 
is represented by the diffusion current ADq / Ld- It must 
be equal to the current on the left, which is the sum of 
the diffusion and the source terms —AD{T)/Ld — S{T). 
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FIG. 2: Illustration of the setup for spin current generation. 
The front gate is close to 2DEG. This gate controls the spin- 
orbit coupling constant via periodic oscillations of Vi. 4> de- 
notes the phase shift between electric potentials of the front 
and back gates. 



The factor A is thus 



A 



p{T)T 



da 



2m*ao [Do + D{T)) dT 



(45) 



Using the expression D = VpT/2 and p = kp/2TT we 
obtain the polarization at the right edge 



py = _ 



m* p{T) 



da 



Trhaa{p{T) + po) dT' 



(46) 



We restored a conventional dimensionaUty by writing h 
in a proper place. An important feature of this expres- 
sion is that py does not depend on the absolute value 
of the SOI coupling constant, but rather from the rela- 
tive amplitude of its variations in time. Hence, the same 
spin injection effect is expected for both, QW with the 
large Rashba SOI, and QW with not large SOI, provid- 
ing that the relative variations of a in time are the same 
in both cases. Another interesting phenomenon, which 
follows from 1)46(1 . is that the AC modulated SOI can re- 
sult in stationary spin accumulation near the gate edges. 
This is due to time dependence of the electron density 
in Eq. H46|l . The DC spin density is obtained by time 
averaging of this expression. A phase shift between p(T) 
and a{T) is required to have this average nonzero. Such 
a shift can be achieved, for example, by manipulating 
the front and back gates, as shown in Fig. |21 Both the 
back and the front gates are equally efficient to induce 
electron density oscillations. At the same time, as shown 
by Grundler,^ the front gate close to 2DEG is necessary 
to control efficiently a. Therefore, the required phase 
shift can be obtained by choosing appropriate phases of 
Vi and ¥2- 



2. w<s:Ld 

Another interesting regime is the case of narrow gate 
with size w (in x direction) much shorter then the charac- 
teristic diffusion length Ld. Outside the gate region the 



solution of Eq. H43|) has the same form as in the previous 
example, namely, 



(47) 



py = zb^exp --^ 
Ld 



where the - and -I- signs refer to x > and a; < regions, 
respectively. Inside the gated region, as seen from Fig. 2 
near x = 0, the polarization varies very fast. Therefore, 
within this range one can retain in Eq. H43() only terms 
containing derivatives of P^ and 5*. This gives 



-Dix) 



Qpy 
dx 



s{:x) = c. 



(48) 



This equation has a form of the spin current conserva- 
tion low. The constant C is, obviously, equal to the 
current just outside the gated region, where S{x) — 0, 
but still X ^ Ld- Hence, from Eq. 1(47(1 one obtains 
C = DqA/Ld- Further, integrating Eq. 1(48(1 we arrive to 
the expression for py in the near vicinity of a; = 



py ^ ^ 

Ld 



dx 



, S{x') 
D{x') 



A 
Ld 



-— dx' 



Do 
D(x') 



(49) 



This expression must coincide with Eq. 147() at w < x <ti 
Ld, namely, with Aeicp{^x/ Ld) ~ A — Ax/ Ld- Substi- 
tuting the latter instead of py into Eq. I(49() and setting 
X ^ 00, one obtains 



A 



1 



1 

Ld jo 
S{x) 



dx 



Dn 



D{x) 



1 



dx 



D{x) 



(50) 



One can neglect the integral in the square brackets 
because it is small by the parameter w/Ld, except of 
special cases when Dq/D{x) becomes very large in some 
points. Expressing D[x) via p{x), as it has been done in 
the previous problem, we arrive to the simple formula 



A = 



1*2 



dx 



da 
df' 



(51) 



Since the integral in Eq. ifET)) is of the order of wda/dT, 
the polarization injected in the case of a narrow gate is 
smaller by a parameter w/Ld than in the previous case 
of the wide gate. That is because of a strong counterflow 
of the diffusion current reducing the effect of the current 
— S pumped by the gate. Also, unlike the previous ex- 
ample, the narrow gate can not inject the stationary spin 
polarization. The rectification effect is absent in this case 
because in Eq. ((51(1 there are no time dependent param- 
eters beside da/dT- 

Figureiniillustrates the spin current [Fig.lSJa)] and spin 
polarization [Fig. Elb)] induced by the narrow AC gate. 
The diffusion current exponentially decays far from the 
gate region with spin accumulation shown in Fig. |2tb) . 
The spin polarization has opposite signs on two sides of 
the gate. 
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(a) 





Assuming that the frequency of a variations is small 
compare to the spin relaxation rate Ts, we obtain the 
general solution of Eq. (|53|l in the form 



(54) 



where k ~ l/Lo- The solution satisfying the boundary 
condition is easily obtained as 



TO 



sinh {ny) da 



2hTTa cosh I 



dT' 



(55) 



For an extended 2D electron system with d — )■ oo, it is 
easy to see that P^{y,T) indicating the absence of 
the bulk spin density, in agreement with the result of sec- 
tion IIIA. At the same time, as expected, the polarization 
is accumulated near y = ±d/2, decreasing exponentially 
when the distance from the boundary increases. 



FIG. 3: Illustration of the spin current flow induced by a 
narrow AC gate: (a) spin current as a function of x, which 
is exponentially decays outside the gate; (b) spin polarization 
as a function of x, which is accumulated around the AC- 
biased narrow gate. Polarization has opposite signs on two 
sides of the gate. 



B. Spin flow in a channel 

We consider a 2D channel (strip) of the width d, so 
that particles are confined in the y direction and free to 
move in the x-direction. The time dependent gate atop 
of this channel causes time variations of the Rashba SOI. 
We assume d ^ / and apply the drift-diffusion equation 
to describe the spin transport in the channel. Since the 
gate is supposed to be infinite in the x direction, the spin 
polarization does not depend on x and only derivatives 
with respect to y and T have to be retained in Eq. (|40|l . 

In this case one can look for a solution of Eq. H46|) 
such that py = P"^ = and = P'^iy). It is easily to 
check that with such choice of P the spin current Eq. 142|) 
in the y direction docs not contain y and z polarization 
components. The homogeneous y-polarized current given 
by the last term in Eq. (|40() fiows only in the x direction. 
Hence, the only current fiowing along the y axis is x- 
polarized. For this current the boundary condition is 
imposed that {J'y)x=±d/2 = 0- Since the time dependent 
a in the last term of Eq. (|42() generates the spin current 
Jy — {l/2)pTda/dT, to satisfy the boundary condition 
this current must can be compensated by counterfiow 
of the spin diffusion current given by the first term in 
Eq. H42fl . Hence, the boundary condition is 



y dT 



0. 



x = ±i 



In its turn, the equation for has the form 



Qpx 

'df 



D{T) 



Q2px 



r.P^(2/). 



(52) 



(53) 



V. SUMMARY AND DISCUSSION 

We have employed the electric gate effect on the 
Rashba spin-orbit interaction in narrow gap semiconduc- 
tor QW and considered spin transport in a 2DEG with 
the space-time dependent Rashba spin-orbit interaction. 
The variations of SOI in time and space were assumed 
to be provided by gates of various shapes. Spin trans- 
port was considered in the framework of the Keldysh for- 
malism which has been applied to derive the Boltzmann 
equation for the spin distribution function in the phase 
space. This equation was further employed in derivation 
of the drift diffusion equation for the spin density. We 
found out that besides the usual terms, both the Boltz- 
mann and drift diffusion equations contain the spin cur- 
rent motive force proportional to the time derivative of 
the SO coupling parameter. 

We have considered several examples with various gate 
geometries. Although these examples do not embrace 
many other interesting possibilities, nevertheless, they 
demonstrate how spin current flow and spin density accu- 
mulation can be electrically controlled by means of gates. 
It has been shown that in some geometries the AC bias 
applied to the gate can result in the DC spin current 
and the stationary spin accumulation. This is a simplest 
self-rectification effect. Probably, more efficient could be 
a special rectification setup consisting of several gates 
in series and combinations of back and forward gates to 
control separately electron density and SO parameter. 
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